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Abstract. We present the X-ray spectrum of TWHya observed at high and intermediate spectral resolution with the Reflection 
Grating Spectrometer (RGS) and the European Photon Imaging Camera (EPIC) onboard the XMM-Newton satellite. TWHya 
is the first classical T Tauri star for which simultaneous X-ray data with both high spectral resolution and high sensitivity were 
obtained, thus allowing to probe the X-ray emission properties of stars in the early pre-main sequence phase. Despite TWHya's 
high X-ray luminosity in excess of lO'^" erg/s its X-ray spectrum is dominated by emission lines from rather cool plasma (T 
~ 3 MK), and only little emission measure is present at high temperatures (T ~ 10 MK). We determine photon fluxes for the 
emission lines in the high resolution spectrum, confirming the earlier result from Chandra that the predominant emission is 
from neon and oxygen, with comparatively weak iron lines. Further, the line ratios of He-like triplets of nitrogen, oxygen and 
neon require densities of rie ~ 10^'^ cm~^, about two orders of magnitude higher than for any other star observed so far at 
high spectral resolution. Finally, we find that nearly all metals are underabundant with respect to solar abundances, while the 
abundances of nitrogen and neon are enhanced. The high plasma density, the (comparatively) low temperature, and peculiar 
chemical abundances in the X-ray emitting region on TW Hya are untypical for stellar coronae. An alternative X-ray production 
mechanism is therefore called for and a natural explanation is an accretion column depleted of grain forming elements. The 
metal depletion could be either due to the original molecular cloud that formed TWHya or due to a settling of dust in the 
circumstellar disk of TW Hya. 

Key words. X-rays: stars - stars: individual: TWHya - stars: pre-main sequence, coronae, activity - accretion 

1 . Introduction of discovering wTTS which can easily escape other identifi- 

cation methods because of their lack of conspicuous features 



T Tauri stars (TTS) are late-type pre-main sequence (PMS) ^he optical or near-infrared. Extensive X-ray studies of a 

; stars, usually found in or near star formation regions. The two y^-jg^y ^tar forming regions were undertaken with ROSAT. 

subclasses of TTS, termed 'classical' and 'weak-line', are dis- i„ particular the spatial completeness of the ROSAT All-Sky 

tinguished by the strength of their emission. This empiri- ^^^^^ ^as resulte d in a large number o f newly identified 

cal distinction most likely also reflects different evolutionary PMS stars (e.g. iNeuhauser et al. 19951 lAlcalaetal. 19951 

states: Most of the 'classical' TTS (cTTS), i.e., those with iwichmann et al 1996b 
strong Hq emission, show signatures of circumstellar matter 

(e.g. IR and UV excesses) accreting onto the central star, while jhe origin of the X-ray emission from wTTS is commonly 
in contrast, 'weak-line' TTS (wTTS) show no signs for the interpreted in terms of scaled-up solar-type magnetic activity, 
presence of disks which presumably have been lost. A detailed ^ut the situation is less clear for their accreting counterparts, 
description of early stellar evolution and the observable stel- ^he cTTS. While nothing should prevent the formation of mag- 
lar properties related wifli the di fferent phases can be found in ^^^^^^ flej^s and ensuing coronae in cTTS, X-rays could alter- 
[Feigelson & Montmerle (1999)| natively also be produced above the hot spot where a magnet- 
X-ray emission from TT S was first detected with ically funnelled accretion flow impacts on the surface of the 
the Einstein satell i te (e.g. [Feigelson & DeCampli 1981| star JLamzin 1999^ . A distinction between a corona or an ac- 
IWalter & Kuhi 19811 IMontmerl e et al. 1983> . It was soon cretion funnel as the site of X-ray emission in cTTS requires a 
recognized that X-ray observations provide an efficient means precise assessment of the physical conditions in the X-ray emit- 
ting region. With the lower resolution X-ray data available from 
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though some differences between the two types of TTS were 
indicated by statistical comparisons of their emission level and 
variability (see e.g. Neuhai lser et al. 19^ IStelzer et al. 20001 
IStelzer & Neuhauser 200 lal . 

With the advent of high-resolution X-ray spectroscopy 
with the grating spectrometers onboard the XMM-Newton 
and Chandra observatories this situation has fundamentally 
changed. The reflection gratings onboard XMM-Newton and 
the transmission gratings onboard Chandra have spectral reso- 
lutions sufficiently large to resolve individual emission lines in 
the soft X-ray range. These emission lines and in particular the 
analysis of their relative strengths provide unprecedented tools 
for the plasma diagnostics of the outer envelopes of stars and 
allow an assessment of the temperature structure, density, and 
chemical abundances in the X-ray emitting region. 

The high signal-to-noise ratio (SNR) required for a mean- 
ingful analysis of high-resolution X-ray spectra clearly fa- 
vors the study of the nearest star forming regions. The 
nearby TTS TWHya (located at a distance of only ~ 57pc; 
IWichmann etaTTg gS I was first thought to be an isolated cTTS 
( IRucinski & Krautter 1 983 1, but later five more stars were 
found in the same region and shown to be members of a physi- 
cal cluster of young stars dominated by TW Hya. The whole 
group is now known as the TW Hydra association (TWA; 
IKastner et al. \991i . which is currently known to host about 
20 young stars (youth being probed by Li absorption), most 
of which - in contrast to TW Hya itself - are wTTS. 

The first - and so far only - example of a high-resolution X- 
ray spectrum of an accreting PMS star was a C/zant/ra/HETGS 
spectrum of TWHya presented by ' Kastner et al. (2002)1 
Interestingly, flux ratios of the OVII triplett lines indicated den- 
sities much higher than found in any stellar corona before, 
and consequently Kast ner et al. (2002)| called into question the 
solar-stellar analogy for TW Hya and cTTS in general. About 
a year after the Chandra observation TWHya was observed 
by XMM-Newton, resulting in a higher SNR grating spectrum 
with the XMM-Newton Reflection Grating Spectrometer (RGS) 
as well as a simultaneously measured medium-resolution CCD 
spectrum with high sensitivity over a broad energy range {E — 
0.3 — lOkeV) from the European Photon Imaging Camera 
(EPIC). The purpose of this paper is to present the results of this 
30ksec XMM-Newton observation and properly discuss their 
implications in the context of X-ray emission from TTS. 

2. Observations 

TW Hya was observed with XMM-Newton on July 9, 2001 for 
a duration of 30ksec (obs-id 01 12880201) with the RGS as 
prime instrument. The observations were performed in full- 
frame mode employing the medium filter for both the MOS 
and pn cameras of EPIC. Because of the optical brightness 
of TWHya (V = 11 mag) the Optical Monitor (OM) was in 
blocked position so that no concurrent optical data is available. 
The observing log for all X-ray instruments onboard XMM- 
Newton is given in Table[2 

TW Hya is by far the brightest X-ray source and the only 
known member of the TW Hya association in this EPIC field. 
By inspecting the time-resolved background emission we ver- 



Table 1. Observing log for the XMM-Newton observation of 
TWHya on July 9, 2001. "Start" and "Stop" refer to start and 
end of exposures. 



Instr. 


UT [hh:mm:ss] 


JD - 2452099 


Expo 




Start 


Stop 


Start 


Stop 


[ksec] 


pn 


06:35:38 


13:57:31 


0.774306 


1.081250 


26.52 


MOS 


05:58:02 


14:01:48 


0.748611 


1.084028 


28.98 


RGS 


05:51:39 


14:05:15 


0.743750 


1.086806 


29.64 



ified that the background was low and showed only little time 
variability without any significant solar-activity related con- 
tamination, so that the full observing time provided data use- 
ful for analysis. The pn data shows slight evidence for pile-up. 
Therefore, we base our analysis on the MOS detectors, which 
furthermore provide a better energy resolution. 

We constructed an X-ray light curve of TW Hya by bin- 
ning all photons from a circular area around the target position. 
We used an extraction radius of 40" including about 90 % of 
the source photons. For the background subtraction we used 
a source-free circle of the same area located near TWHya. 
Background subtracted MOS lightcurves of TW Hya in the en- 
ergy range 0.3— 2keV are shown in Fig.^ The counts collected 
by MOS 1 and MOS 2 have been added to increase the SNR. 

To investigate spectral changes during the observation we 
defined two energy bands, S for 'soft' and H for 'hard', and 
computed a hardness ratio HR = {H — S)/{H + S). The 
overall X-ray spectrum of TW Hya is very soft, with only a 
few photons recorded above ^ 2keV (see discussion below). 
To ensure sufficient signal in both bands we chose S = 0.3 — 
1.0 keV and H = 1.0 - 2keV. The time evolution of HR is 
shown in the lower panel of Fig.^ Neither the Ughtcurve nor 
the evolution of the HR show conspicuous signs for flaring. 



3. High-resolution Spectrum: RGS 

The RGS data of TW Hya were analysed using the standard 
pipeline rgsproc of the XMM-Newton Science Analysis System 
(SAS), version 5.4.1. We extracted the total (source plus back- 
ground) and a background spectrum for TW Hya. In this pa- 
per only the first dispersion order is discussed. The full first 
order RGS spectrum of TWHya is displayed in Fig. |2l ex- 
cept for the small range between 35 — 37 A where no spectral 
lines were found. In Fig.|2lidentifications of the stronger emis- 
sion lines are also provided from comparisons with line lists 
of the MEKAL code JMewe et al. l"985llMewe et al. 1995> . the 
National Institute of Standards (NIST) data base, and the re- 
sults of solar flare observations (Phillips et al. 1999) . In the 
RGS spectrum one recognizes the lines of H-Uke and He-Uke 
carbon, nitrogen, oxygen and neon, while iron lines are weak 
or absent (especially those near 15 and 17 A). The typical sil- 
icon and magnesium lines below 10 A are also absent in the 
RGS spectrum. Note, however, that at those wavelengths the 
RGS effective area decreases rapidly. 
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Fig. 1. Background subtracted MOS lightcurves of TW Hya in 
the energy range 0.3 — l.OkeV and 1.0 — 2.0keV. The panel 
on the bottom shows the hardness ratio computed from these 
bands; binsize is 1 ksec. 



3. 1. Identification and Analysis of Emission Lines 

To determine fluxes for individual emission lines we used the 
CORA' line fitting tool. CORA is based on a maximum likeli- 
hood technique (see N ess & Wichmann 2002> and fits a spec- 
ified number of emission lines combined with a vertical offset 
to account for a constant background (both instrumental and 
due to continuum emission from the source) assuming some 
specified emission line profile function. Earlier studies have 
shown that a Lorentzian profile provides a reasonable descrip- 
tion of the shape of RGS line profiles (e.g. Ness et al. 20021 
IStelzer et al. 2002 1. Before feeding the spectra into CORA we 
corrected for the different extraction area of total and back- 
ground spectrum. To take into account possible small back- 
ground variations we selected successive wavelength intervals 
in which the background can be approximated by a constant 
and performed the analysis separately for each of these inter- 
vals. Some ambiguity does remain in the determination of the 
local background level. To account for uncertainties associated 
with the background level we add a systematic error to the sta- 
tistical error of the line counts. This systematic error was de- 
rived by using the highest and lowest background levels that 
would provide an acceptable line fit after visual inspection of 

' CORA can be downloaded from 

|http://www.hs.uni-hamburg.de/DE/Ins/Per/Ness/Cora 



the solution. Finally, the line counts and their errors (both sta- 
tistical and systematic) were transformed into photon fluxes 
with the help of the effective areas of the RGS. 

A list of the brightest lines in the spectral range of the RGS 
(6 — 37 A) is provided in Tabled where we list the lines' cen- 
tral wavelengths and their errors, the fitted number of counts 
in the line including errors and the derived photon fluxes. A 
suggested identification is also given. A substantial part of the 
total wavelength range is covered only by one instrument due 
to failure of one CCD on each RGS, and because of the gaps 
between adjacent CCDs. For lines in these regions we give the 
result of the active RGS, for the lines covered by both RGS we 
summed up the two spectra. For diagnostic purposes the most 
important lines are the Lya and triplet lines of hydrogen-like 
and helium-like ions. The measurements of these lines for all 
species detected in the RGS spectrum and in addition the most 
prominent iron lines are provided in Table |2l We emphasize 
again the unusual weakness of iron lines as one of the most 
remarkable features of the X-ray spectrum of TW Hya. 



3.2. Line Temperatures 

For the elements neon, oxygen and nitrogen both the Lya 
emission lines at 12.1 A, 19.0 A, and 24.8 A respectively, and 
the He-like resonance lines at 13.5 A, 21. 6 A, and 28.8 A re- 
spectively are detected in our RGS spectra. In Fig.|3]we plot 
the expected energy flux ratio from the H-like Lya and He- 
like resonance lines for neon (dotted line), oxygen (solid line) 
and nitrogen (dashed line) versus logarithmic temperature us- 
ing the CHIANTI calculations dDere et al. 19971 . We also in- 
dicate the measured line ratio range for neon (upright shad- 
ing), oxygen (right tilted shading) and nitrogen (left tilted shad- 
ing). These line ratios clearly correspond to effective "isother- 
mal" temperatures, which, however, are independent of abun- 
dance. For oxygen and nitrogen these temperatures agree at 
IgT [K] « 6.45, for neon a slightly larger effective tempera- 
ture is found. The X-ray spectrum of TW Hya does not contain 
any strong high temperature lines despite its large overall X-ray 
luminosity. Thus we conclude that most of the observed line 
emission seems to be produced by material located at pretty 
much the same temperature. In fact, we verified that a 1-T ther- 
mal model (VMEKAL) with free abundances can provide an 
acceptable description of the RGS spectrum (Xrcd ~ ^-^^ 
2333 d.o.f.). 

However, one needs to keep in mind that the RGS effective 
area decreases rapidly at wavelengths below ~ 10 A, i.e. en- 
ergies above 1.2keV. Elements such as Mg and Si have their 
characteristic He-like and H-like line transitions in this wave- 
length region. In the RGS spectrum of TW Hya these lines re- 
main undetected because of lacking SNR, but this critical spec- 
tral region can be studied with EPIC at medium resolution (see 
Sect.l^. The EPIC spectrum does also show a high-energy tail 
that can not be explained by the rather low plasma temperatures 
dominating the RGS line emission (cf.. Sect. Sand Fig. 13. 
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Fig. 2. First order background subtracted XMM-Newton RGS count rate spectrum of TWHya. Overplotted is a 3-T VMEKAL 
model whose parameters were derived from a combination of emission line analysis and global fitting of the medium-resolution 
MOS spectrum. The MOS spectrum is shown together with the same model in Fig.|3 and the best fit parameters are summarized 
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Table 2. Spectral lines identified in the first order XMM-Newton RGS spectrum of TW Hya. left - Observed line center, width, 
number of counts, and photon flux. Uncertainties include statistical 1 a error plus a systematic error stemming from the uncertainty 
of the local background level (see text for details), right - Line identifications from |Mewe et al. (1985) Mew e et al. (1995)| and 
[Phillips et al. (1999)| For line blends we list the most likely contributors according to X-ray observations of other stars. 



A 


a 


J 


Photon Flux 


Identification 




[A] 


[A] 


[cts] 




Ion 


Trans. 


X 


RGSl 


18.983 


0.055 


465.7 ±43.7 


30.4 ± 2.8 


OVIII 




18.97 


21.616 


0.047 


241.6 ±28.4 


17.4 ± 2.1 


OVII 


r 


21.61 


21.836 


0.053 


128.2 ± 20.1 


9.5 ±1.5 


OVII 


i 


21.80 


22.100 


0.037 


5.3 ±7.1 


0.4 ±0.6 


OVII 


f 


22.10 


RGS 2 


12.152 


0.060 


111.4 ± 26.1 


6.3 ± 1.5 


NeX 




12.13 


13.460 


0.050 


402.0 ± 38.1 


23.1 ± 2.2 


NelX 


r 


13.45 


13.560 


0.057 


308.9 ± 19.9 


17.4 ± 1.1 


NelX 


i 


13.55 


13.710 


0.054 


149.7 ± 27.3 


8.2 ± 1.5 


NelX 


f 


13.70 


17.065 


0.075 


46.8 ± 15.9 


2.7 ±0.9 


FeXVII 


NeSG 


17.06 


24.793 


0.060 


171.5 ±31.5 


11.8 ± 2.2 


NVII 




24.78 


28.820 


0.059 


32.1 ± 14.2 


2.9 ± 1.3 


NVI 


r 


28.79 


29.101 


0.030 


14.2 ± 7.2 


1.3 ±0.6 


NVI 


i 


29.10 


29.540 


0.030 


4.3 ± 5.8 


0.4 ±0.5 


NVI 


f 


29.54 


RGS 1+ RGS 2 


16.030 


0.056 


118.7 ±37.2 


3.5 ± 1.1 


Fe XVIII 




16.00 










OVIII 




16.01 


18.649 


0.037 


53.4 ± 12.1 


1.7 ±0.4 


OVII 


He3A 


18.63 


33.757 


0.065 


113.9 ± 20.7 


7.8 ± 1.4 


CVI 




33.70 



1 00.00 F 



.9 IO.OOe- 



1 .00 E- 



0.10 r 



0.01 



6.1 




6.40 6.50 

Log T 



6.80 



Fig. 3. Temperature dependence of the Lya/r line ratio for 
neon (dotted line), oxygen (solid line), and nitrogen (dashed 
line) according to the CHIANTI code (Deree tal. T997T >. The 
ratios observed in the RGS spectrum of TW Hya are indicated 
by the shaded patterns. All observed values are compatible with 
a 2 — 3 MK hot isothermal plasma. 



3.3. Differential Emission Measure 

In order to quantify a differential emission measure (DEM) 
distribution we adopted the following approach; see also 
ISchmitt & N ess 2004l for more details on this DEM reconstruc- 
tion method. We first introduced a fixed temperature grid lo- 



cated at grid points log T = [6.3,6.5,7.3]. The two lower tem- 
peratures are based on the result from the above described line 
analysis, while the highest temperature is required for a match 
to the medium-resolution EPIC spectrum whose high-energy 
tail cannot possibly be described with a 3 MK plasma alone. It 
turns out, however, that the contribution of the hot component 
to the overall RGS count rate is very small. 

At each temperature grid point Ti (i=l,3) we assume the 
presence of some emission measure EMi. Denoting the el- 
ements carbon (i=l), nitrogen (i=2), oxygen (j=3) and neon 
(j=4) with the index j, we can then compute the line emissivi- 
ties Ph,!.} and Pncij for the hydrogen- and helium-like lines of 
the elements j for the temperature components i. The observed 
line ratios pobs.j (j=2,4) can then be modelled as 



Pmod.j 



^f^i EM, ■ FhcU 



(1) 



and we seek a best fit set of the emission measures EAIi 
through minimizing 



X 



4 

E 

j=2 



(PobsJ Prnod,] ) 



(2) 



where crj denotes the error of the line ratio j. Clearly, since 
we are dealing with line ratios of the same elements, the nor- 
malization can be set arbitrarily (we chose EMi = 1). The 
thus derived EM-distribution is abundance independent. The 
best fit values for the (relative) values of EMi are found to be 
[1.0,2.35,0.59]. This choice of emission measure reproduces 
the observed line ratios extremely well (x^ed ~ However, 
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Table 3. Best fit parameters for the X-ray spectrum of TW Hya 
described by a multi-temperature thermal model (VMEKAL) 
including photo-absorption (N^ = 2 x 10^" cm^^). The elec- 
tron density was fixed to 5 x 10^^ cm^^. Elements not listed 
below are irrelevant in this spectrum, and their individual abun- 
dances were held fixed at the solar value. The temperatures and 
the coupling between some parameters were pre-determined 
by the analysis of the emission lines seen in the RGS spec- 
trum. Values for the abundances are given with respect to the 
solar photosphere, logT is in K, and log EM in cm~^. (A) 
- Best fit 3-T model to the EPIC MOS spectrum in the range 
A = 5 - 40 A. (B) - Same as (A) with free Mg and Si abun- 
dances. 



-Model 



-Model B- 



Parameter 



par am. 
coupling 



param. 
value 



param. 
coupling 



param. 
value 



C 
N 
O 

Ne 
S 

Mg 

Si 

Fe 



■ 1.07 X O 

■ 3.14 X O 



= 0.28 
= 0.82 
0.26 
3.15 
= 0.0 
= 1.0 
= 1.0 
0.26 



1.07 X O 
■ 3.14 X O 



0.40 
1.18 
0.37 
3.97 
= 0.0 
0.65 
0.13 
0.24 



logTi 
logTa 
logTs 



6.3 
6.5 
7.3 



6.3 

: 6.5 
: 7.3 



log EMi 
log EM2 
log EM3 



2.35 X EMi 
0.59 X EMi 



52.64 
53.01 
52.41 



2.35 X EMi 



52.53 
52.90 
52.53 



Xrod (d.O.f.) 



2.79 (233) 



2.61 (232) 



1.00 r 



0.10 ^ 



0.01 




given the ill-conditioned nature of the reconstruction problem, 
other choices of the EM-distribution also provide statistically 
acceptable fits. With the thus derived set of best fit EMi values 
we can choose the normalization to reproduce, say, the oxy- 
gen Lija fine fluxes. Then we use this normalization to deter- 
mine abundances of carbon, nitrogen, neon and iron relative 
to oxygen by forcing agreement between predicted and mea- 
sured fluxes in the observed lines. Note that for carbon and 
iron we can only use the Ly^ line (for carbon) and the Fe XVII 
17.07 A line (for iron). Using an interstellar hydrogen column 
density of Nh = 2 x 10^" cm~^ - which is required for a good 
fit of the EPIC data and is in agreement with the expected inter- 
stellar absorption towards TW Hya - we find with this proce- 
dure flie values C/0 = 1.07, N/O = 3.15, Ne/O = 9.77 and 
Fe/0 — 0.91 with respect to cosmic abundances. Realistic 
abundance errors are difficult to determine. We estimate abun- 
dance errors of about 20% for neon, and about 30 % for carbon, 
nitrogen and iron. 

4. EPIC medium resolution spectrum: EPIC 

The count rate spectrum of each EPIC MOS detector was ex- 
tracted from the area described in Sect. |2l We extracted and 
analysed the EPIC spectra with SAS version 5.4.1. Individual 



.Fig. 4. EPIC MOS 1 and MOS 2 spectrum of TWHya. The 
data is overlaid by the best fit 3-T VMEKAL model B from 

-Table|3l see also Fig.|2]where the same model is folded with the 
RGS response. The individual contributions of the two lower- 
temperature components and of the high-temperature compo- 
nent are also shown. 



redistribution matrices and an ancilliary response file were gen- 
.erated in the course of our data reduction process. The combi- 
nation of these two files fully models the instrumental response. 

The EPIC spectrum has a pronounced peak at £' 
-0.91 keV or 13.6 A. Since we know from the XMM-Newton 
RGS spectrum that both the continuum and iron line emis- 
sions are weak, we conclude that this spectral bump must 
"represent emission from H- and He-like neon. This re- 
"markable feature is not observed in EPIC data of other 
young stars (see e.g. Giide l et al. 200 1 1 IStelzer et al. 20021 
IStelzer & Neuhauser 2001b 1. emphasizing the peculiar nature 
of the X-ray spectrum of TW Hya. 

Approximating the MOS spectra by a sum of three ther- 
mal models with fixed values of temperature and based on the 
parameters derived from the line analysis, we verify that the 
temperature and abundance structure derived from the analy- 
sis of individual emission lines also provides a consistent de- 
scription of the medium-resolution data. In this global fit to the 
MOS spectrum an absorption term with atomic cross-sections 
of .Morrison & McCammon (1983). and TVh = 2 x 10^° cm-^ 
is applied to describe the effects of inter- and/or circumstel- 
lar cold gas. The electron density was set to 5 x 10^^ cm~'^ 
following the result of our line-based density analysis (see 
Sect. 15. 2^ and kept fixed. We restricted our modeling to the 
5.5-40A range, and used the VMEKAL code as implemented 
in XSPEC, version 11.2.0. 

We first allowed only the oxygen abundance and the abso- 
lute value of the emission measure to vary, i.e., the parameters 
not constrained by our line modeling. The best fit determined 



in this fashion results in Xrcd ~ •^■'^ (^35 d.o.f.), and the resid- 
uals suggest that a slight adjustment of the abundances may be 
needed. If we allow the neon and iron abundances to de-couple 
from the oxygen abundance, the fit improves to Xrcd = 2. 8 
(233 d.o.f.). Note that these "new" abundances are still within 
the errors of the abundances derived from our hne-based RGS 
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analysis. The best fit parameters of this model are summarized 
in Table|3l(Model A). We point out that we do not derive nitro- 
gen and carbon abundances from the EPIC spectra, because the 
instrument has little sensitivity in the relevant spectral range. 
These abundances are solely derived from the RGS data and - 
for carbon - are based only on one single line. 

Up to this point our fit results are in good agreement 
with the findings presented by |Kastner et al. (2002)| from their 
analysis of the Chandra HETG spectrum of TWHya, i.e., a 
strong over-abundance of Ne and emission dominated by low- 
temperature plasma. Inspection of the EPIC fit residuals does 
however show clear discrepancies between the EPIC data and 
model at A < 14 A, where the H-like Lya and He-like triplet 
transitions of magnesium and silicon are located. In particular, 
two bumps in the model come from emission of Si XIV and 
Mg XII without any counterpart in the data. Therefore we al- 
lowed for adjustment of the Mg and Si abundances in the spec- 
tral fitting process. Since any changes in the Mg and Si abun- 
dances will be compensated for by changes in EM3, we must 
also let vary EM3. The resulting best fit models have magne- 
sium and silicon abundances well below solar values with a 
best-fit Xrcd °f 2.61. The 90 % error contour of the magnesium 
is at 0.94 (w.r.t. solar), and the silicon abundance is definitely 
subsolar. A summary of all parameters of this model is given in 
Table |3l (Model B). In Fig.|4|we show the best fit 3-T model 
together with the MOS data. For illustration the sum of the 
two low-energy components and the high-energy component 
are also plotted separately. 

In order to point out the effects of changing Si and Mg 
abundances we consider only the range between 5.5 — 12 A. In 
Fig. Is] we confront the two models from Table |3]in this wave- 
length interval with each other. Clearly the fit with free abun- 
dances describes the EPIC data far better than assuming solar 
abundances for silicon and magnesium. In the 5.5 - 12 A m- 
terval the Xrcd decreases from 1.23 (102 d.o.f.) for the model 
with solar Mg and Si to 0.90 (101 d.o.f.) for the model where 
Mg and Si are free parameters, leading to subsolar abundances 
for both elements. 

As a final consistency check of our procedure we now 
folded back the best fit model B with the response of the RGS. 
The result is overplotted on the RGS data in Fig. |2l As ex- 
pected, we find good agreement between model and data. The 
emission measure of the high-T component is about a factor of 
three lower than the sum of the two low-T components. This 
together with its high temperature leads to a very low contribu- 
tion to the RGS count rate of 0.015 cps (RGS 1) and 0.021 cps 
(RGS 2), implying that the high-T component is essentially un- 
detectable in the RGS. (The lower count rate in RGS 1 is caused 
by the missing CCD in the 10 — 14 A interval.) We thus con- 
clude that the model B described in Table |3] provides an ac- 
ceptable description of both the XMM-Newton EPIC and RGS 
spectra. 

5. Discussion 

So far TWHya is the only cTTS, i.e. a PMS star with strong 
Hq emission, for which a high-resolution XMM-Newton X-ray 
spectrum is available. Since in young stars H^ emission is re- 
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Fig. 5. Close-up of the EPIC MOS 1 and MOS 2 spectrum of 
TW Hya in the range 5.5 - 12 A: fo/^ - Model (A) from TablegJ 
Mg and Si abundance fixed to solar; bottom - Model (B) from 
Table|3l Mg and Si abundance free parameters resulting in sub- 
solar abundances for both elements. 



lated to accretion TWHya presents a unique test case for X-ray 
emission in this early phase of stellar evolution. Furthermore, 
it is one of the few stars observed with both XMM-Newton and 
Chandra at high spectral resolution. In the following we dis- 
cuss the properties of the X-ray emission of TW Hya, and we 
present a physical scenario explaining it. 

5.1. Variability 

During the Chandra observation of TWHya discussed by 
[Kastner et al. (2002)] a sharp rise in intensity followed by a 
nearly linear decay occurred, that the authors interpreted as a 
flare. We find little evidence for variability in the XMM-Newton 
observation presented here (acquired ~ 1 yr after the Chandra 
data). In any case the lightcurve does not carry the clear signa- 
tures of a flare, and no change in spectral hardness is observed. 
Similarly [Kastner et al. (2002)| found no changes in the spec- 
tral parameters during the variable phase of the Chandra ob- 
servation, very untypical for stellar (flare) activity. The flare- 
like event in the Chandra observation has a decay time of 
~ lOksec. Assuming a flare loop at T '--^ 3MK (the dominat- 
ing temperature derived from the X-ray spectrum) that decays 
through radiative cooling the density should be ^ lO^*' cm^"^. 
This is in sharp contrast to the result from the analysis of 
line ratios (see Sect. 15. 2> . according to which the densities are 
higher by two orders of magnitude. 

We directly compare the emission levels by computing X- 
ray fluxes for the time-averaged spectrum. From the Chandra 
HETGS spectrum [Kastner et al. (2002)| derived an X-ray flux 
of 2.5 X 10~^ photons cm"^ s"^ for the 0.45 - 6.0 keV band. 
Integrating the RGS spectrum in the energy band which comes 
nearest to the spectral range of the HETGS, i.e. 5.5 — 27.5 A 
or 0.45 — 2.25 keV, we arrive at a similar value of /x,rgs i ~ 
/x,RGS2 = 2.3 X 10^'^ photons cm^^s^^ corresponding to 
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Fig. 6. Comparison of photon fluxes for individual emis- 
sion lines in the X-ray spectrum of TW Hya measured with 
the XMM-Newton RGS (see Table |3 and the Chandra 
HETGS-MEG (see IKastner et aLIHUa . The XMM-Newton 
and Chandra observations are separated by about one year. All 
lines are constant within about a factor of two. 



-^x.RGS = 1.4 X 10^" erg/s at a distance of 57 pc. This suggests 
that TW Hya exhibits only moderate long-term variability. 

Based on the small amplitude of long-term variations dis- 
cussed above, comparing fluxes measured for individual emis- 
sion lines allows to check the calibration of the instruments 
on both satellites. We observe a trend towards slightly higher 
fluxes in the RGS as compared to the HETGS (Fig. |6}. 
However, the difference is a factor two at most. 

5.2. Density in tine X-ray Emitting Region 

The high-resolution X-ray spectrum obtained with the RGS 
allows us to measure flux ratios between individual emission 
lines and use them for plasma diagnostics. The triplets of He- 
like ions are of particular importance, because ratios between 
the intensity of resonance r, intercombination i, and forbid- 
den / line are sensitive to (electron) temperature (through 
G = (/ + i)/r) and density (through R = J /i). While the G- 
ratio is of Umited value for temperature diagnostics, the i?-ratio 
constitutes a unique tool to probe the electron density in stellar 
coronae. Although the triplets from silicon (Si XIII) to nitro- 
gen (NVI) fall in the spectral range of the RGS, in TWHya 
only three elements - nitrogen, oxygen, and neon - show suffi- 
cient emission in the respective transitions of their He-like ions 
to warrant analysis. 

In Fig.0we zoom in on the spectral region around the oxy- 
gen and neon triplets of TW Hya overlaying the data with our 



line fit result. For the modeling of both triplets we used three 
Lorentzians representing the r, i, and / lines. Some caution 
is in order when interpreting the spectrum near the position of 
the NelX triplet: the triplet is not fully resolved in the RGS, 
and is known to form a blend with a number of iron lines, pre- 
dominantly a Fe XIX line at 13.53 A that contaminates the Ne 
intercombination line. However, we performed fits of the neon 
triplet including iron and find undistinguishable results. Since 
we find little hints for the presence of Fe XIX (and generally 
any iron) in the RGS spectrum (e.g. the lines at 13.80A and 
14.66 A are not seen) we assume that iron contamination plays 
only a minor role and can be neglected. 

The remarkable feature of the He-like triplets of TW Hya 
is the weakness of the forbidden lines. In the LETGS and RGS 
spectra of all late-type stars observed so far the forbidden lines 
are always stronger than the intercombination line (in oxygen 
and neon) and in many cases the measurements are consis- 
tent with the so-called low-density limit (e.g. lNess et al. 200T1 
INess et al. 20021 IStelzer et al. 2002 ». This is definitely not the 
case for TW Hya. We note in particular that with TW Hya be- 
ing rather cool (spectral type K7; iRuc inski & Krautter 1983> 
the low intensity of the forbidden line cannot be attributed to 
UV radiation at the wavelengths corresponding to the / — > i 
transitions JCosta et al. 2000l l. Instead it is a direct indicator 
of high densities. The formally derived //i-value for OVII 
is i?ovii < 0.16, i.e., only an upper limit can be given due 
to the weakness of the OVII /-line. Comparing our measure- 
ments to recent calculations by ' Porquet et al. (20()T)] without 
any radiation field we find a lower limit on the density of 
ne,ovu > 1 X 10^^ cm^'^. This is corroborated by the low flux 
of the nitrogen forbidden line which is only marginally detected 
(cf., Tableill. The i?-ratio for neon (i?Ncix = 0.47 ± 0.11) in- 
dicates a density of rie, NcIX ~ 1 x 10^'^cm^''. Inclusion of 
the Fe XIX blend in the fit of the neon triplet would lower the 
strength of the i line, and lead to a somewhat higher value for 
-RncIX (0.62±0.12), and subsequently lower density. However 
the numbers are compatible with the ones cited above within 
the statistical uncertainties. 

The densities derived from the oxygen and neon triplets are 
at least two orders of magnitude above typical coronal densi- 
ties. Our density measurements thus fully confirm the results 
derived by [Kastner et al. (2002)| from their HETGS spectrum, 
and we conclude in particular that the high densities cannot be 
a result of the flaring reported by iKastner et al. (2002)| in their 
Chandra data. 



5.3. Elemental abundances 

The XMM-Newton RGS spectrum as well as our spectral 
fits to the EPIC spectrum show the X-ray emission of 
TWHya clearly dominated by neon emission in agreement 
with [Kastner et al. (2002)] Our results specifically suggest a 
neon abundance significantly enhanced with respect to the so- 
lar value (A^e ~ 4A^e0), while iron and oxygen, the ele- 
ments with the largest number of lines in the sensitive re- 
gion of the XMM-Newton and Chandra instruments, are un- 
derabundant with respect to solar values, again in agreement 
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Fig. 7. He-like triplets of Ne IX and O VII measured with the RGS together with best fit Lorentzians. The continuum close to 
each triplet is approximated by a straight horizontal line. The intensity enhancement to the right of the Ne IX triplet in the RGS 
spectrum is due to FeXVII. But the major contamination of the NelX triplet is expected to come from a FeXIX line which is 
unresolved from the triplet intercombination line (see text in Sect. l5.2l l. 



with [Kastner et al. (2002) For the elements silicon, magne- 
sium, anTliittogeir|Kastne^ral. (2002) claim solar abun- 
dances, while carbon lines were not covered by the HETGS 
spectrum. Interestingly, the lines of He-like and H-like silicon 
are clearly visible in the HETGS spectrum of TWHya (cf.. 
Fig. 2 in lKastner et al. 20028 . while the same authors note "a 
curious lack of Mg features, relative to the [i.e. their] model." 
While the RGS spectra are not sufficiently sensitive in the rel- 
evant spectral region, the EPIC spectrum clearly suggests sub- 
solar abundances also for magnesium and silicon. 

Let us now focus on the line emission from neon, oxy- 
gen, nitrogen, and carbon as seen by the RGS. Our spectral 
analysis both with RGS and EPIC shows that very little low 
temperature material (with logT [K] < 6.3) can be present. 
|Ness et al. (2002)| report photon flux ratios for Lya/r of 0.59 
(oxygen) and 0.94 (nitrogen) for the low-temperature F-star 
Procyon, and 1.87 (oxygen) and 3.26 (nitrogen) for the higher- 
temperature K-star e Eri, which must be compared to the values 
of 1.75 (oxygen) and 4.07 (nitrogen) measured for TW Hya de- 
spite its overall low temperature. Although DEM distributions 
are typically not unique because of the ill-defined nature of 
the inversion problem, we want to emphasize that the peculiar 
chemical abundances of TW Hya are not caused by our choice 
of the EM distribution. For this we examined the energy flux 
ratio between the Ly^ lines of C VI and N VII as a function of 
temperature and find that this line ratio can not be explained by 
any single temperature or combination of temperatures under 
the assumption of solar abundances for both elements (Fig.|8]a). 
Therefore nitrogen must be overabundant relative to carbon 



(and also oxygen) with respect to the solar abundance ratio, 
again in accordance with our fit results. We then investigated 
the flux ratio between the r line of neon and the oxygen Ly^ 
line in the same way (Fig.|8]b), and again the observed flux ratio 
can not be reproduced by any single temperature or combina- 
tion of temperatures if the abundances were solar, implying that 
neon must be overabundant compared to oxygen with respect 
to solar abundances, as found in our fit results. Finally, inspect- 
ing the flux ratio between the C VI Ly^ line and the O VII r 
line we find perfect consistency assuming the C/O ratio being 
solar. Since our EPIC results suggest that silicon, magnesium, 
iron, and oxygen all have similar abundance (with the possible 
exception of magnesium) we thus conclude that the elemental 
abundances of these elements as well as that of carbon are re- 
duced with respect to solar values by about a factor of 3 — 4. 
Nitrogen is enhanced with respect to these elements by a fac- 
tor of 3 as evidenced by the XMM-Newton RGS spectrum, and 
neon is enhanced by a factor of 10 as evidenced by both the 
XMM-Newton RGS and EPIC spectra. 



5.4. Interpretation 

The three properties distinguishing the X-ray spectrum of 
TWHya from that of "normal" late-type stars are (1) the ab- 
sence of a hot component with large emission measure, (2) the 
high density of the X-ray emitting region, and (3) the peculiar 
elemental composition, and in particular the low abundance of 
iron. The first two properties can be naturally explained by at- 
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Fig. 8. Temperature dependence of various line ratios measured in the RGS spectrum of TW Hya. Solid lines are CHIANTI model 
calculations (Dere et al. 1997i, shaded regions denote the 1 a range observed with the RGS. (a) - N VII Ly^ / C VI Lya\ (b) - 
Ne IX r / O VIII Ly^, (c) - C VI Ly^ I O VII r. The plots demonstrate that the C/O abundance ratio is close to solar, while the 
N/O and Ne/O abundance ratio can not possibly be solar. 



tributing the X-ray emission to an accretion shock, the third by 
the specific environment of TW Hya. 
Employing the strong shock formula 

16fc 

with /i denoting the mean molecular weight, mn the hydrogen 
atom mass, and k Boltzmann's constant, we can convert the 
measured temperature lgr[K] ~ 6.45 to an infall preshock 
velocity V^rc of about 350/ \/ (/i) km/s. The preshock parti- 
cle density nprc should be one fourth of the measured post 
shock density ripost, i-e., 2.5 x 10^^ cm^'^. In order to estimate 
Lpost, the thickness of the postshock cooling layer, we equate 
the infalling kinetic energy flux with the radiated energy flux 
JLamzin 1995> through 

/^"tHnpre^pre _ 2 , j 

2 ~ '^post^'^mcan -impost- 

Here Amoan denotes the total mean radiative post shock cool- 
ing rate per unit emission measure. Denoting Amean in units of 
lO"^'^ erg cm'^ s~^, we derive 

9 X 10^ ,_i . , 

-t^post — ^ — (,Anican,23j Cm . 

For an isothermal cosmic abundance plasma at Ig T [K] = 6.45 
Amcan.23 IS generally close to unity; if we assume A„ican,23 = 
0.3 because of the depleted elemental abundances and /i ~ 1.5 
the post-shock thickness is Lpost ~ 2500 km with admit- 
tedly considerable uncertainty. Using the observed O VIII en- 
ergy flux of 3.2 X 10^^'^ erg/cm^/s and an assumed OVIII 
Lya cooling function of 9.6 x 10^^^ erg cm"^ s^^ (assuming 
lgT[K] — 6.45 and an oxygen depletion of 0.3) we com- 
pute an overall volume emission measure of EM = 1.3 x 
lO^'^ cm^'^ at the distance of 57 pc. Since ripost ~ 10^'^ cm^'^, 
the emitting volume Vomit must be Vomit = 1.3 x 10^^ cm'^. 
The shock area Ashock then becomes ^shock = Vomit /^post = 
5.1 X 10^* cm^, which would constitute only a very tiny 
fraction of the visible stellar surface (< 0.01 %). Clearly, 
these estimates are by necessity very rough, but in quali- 
tative agreement with previous observations in other wave- 
lengths: |MuzerolleetaL_(200^ shown that the shape of 



the optical/UV spectral energy distribution of TWHya can 
be explained by an accretion shock which fills ^ 0.3% of 
the stellar surface, and small accretion shocks are also found 
on other T Tauri stars ( |Calvet & Gullbring 1 998 1. Computing 
the mass accretion rate in the funnel flow through Mace = 
l/4Ashock"-postA^"^//Vpro we find Afacc 1 x 10"" Mo/yrs, 
which is entirely plausible and within one order of magnitude 
of the mass accretion rate of TW Hya derived from UV data 
(jMuzer oll e et al. 200"2t . 

What kind of material is flowing through this funnel onto 
TW Hya ? We have shown that this material is in gen- 
eral metal-depleted in particular with respect to grain-forming 
elements such as Mg, Si, Fe, C and O. Overabundances 
of neon are found in stellar coronae of RS CVn systems 
(lAudard et al. 20031 . and nitrogen overabundances in evolved 
stars JSchmitt & Ne ss 2002> . The low abundances especially of 
iron and carbon call for a different explanation for these chem- 
ical peculiarities. TWHya is rather young and formed "just 
recently" out of a molecular cloud. Models of cloud chem- 
istry at high densities (n > 10^ cm"'^) predict that almost all 
metals condense into grains except nitrogen ( |Charnley 1997^ . 
If the remaining depleted gas is then decoupled from the 
grains in an accretion disk around a young star (as is ex- 
pected in planetary formation models), the accreted and 
shock-heated gas is expected to be metal-deficient as well. 
Indeed, [Weinberger et al. (2002)| found evidence for the pres- 
ence for condensed siUcates in the disk of TWHya, and 
[Herczeg et al. (2002)| speculated about a deficiency of Si in 
the accretion flow of TW Hya based on the examination of its 
far-UV spectrum. XMM-Newton has completed this picture of 
TW Hya as an accreting PMS star, and for the first time pro- 
vided a consistent view of X-ray emission from a stellar accre- 
tion shock by observation. 

6. Summary 

To summarize, the marked differences between the X-ray prop- 
erties of TW Hya and active main-sequence dwarf stars sug- 
gest either (i) strong evolutionary effects on coronal structures, 
or (ii) an entirely different origin for the X-ray emission of 
TW Hya. X-ray emission from a shock at the bottom of an ac- 
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cretion column provides a plausible mechanism for the latter 
possibility. The derived shock temperatures, gas densities, fill- 
ing factors and mass accretion rates as well as the obvious lack 
of grain-forming elements in the X-ray emitting region provide 
strong indications in favor of this scenario. 

A high-resolution X-ray spectrum of another member of 
the TW Hya association, HD 98800, was recently obtained with 
Chandra (CXC Press Release 03-04). The object is a quadru- 
ple star with a circumbinary disk. In terms of Ha emission 
HD 98800 is a wTTS. Therefore it does not come as a sur- 
prise that its X-ray properties are very different from those 
of TW Hya, and clearly reminiscent for typical stellar activ- 
ity. Similarly, the high-resolution X-ray observation of PZ Tel, 
a ^ 20Myr old post-TTS in the Tucanae association, does not 
show properties drammatically different from older late-type 
stars (Argiroffi et al., in prep.). 

We conclude that TW Hya stands out as a presently unique 
example for a star with strong evidence for accretion signa- 
tures in its X-ray spectrum. Future high spectral resolution and 
high sensitivity X-ray observations will have to show whether 
soft emission, high densities, and metal-deficiency are general 
characteristics of cTTS. 
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